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ABSTRACT 

We present narrowband images of the gravitational lens system Q 2237+0305 made with the Nordic Optical 
Telescope in eight different filters covering the wavelength interval 3510-8130 A. Using point-spread function 
photometry fitting we have derived the difference in magnitude versus wavelength between the four images 
of Q 2237+0305. At A = 4110 A, the wavelength range covered by the Stromgren-v filter coincides with the 
position and width of the CIV emission line. This allows us to determine the existence of microlensing in 
the continuum and not in the emission lines for two images of the quasar. Moreover, the brightness of image 
A shows a significant variation with wavelength which can only be explained as consequence of chromatic 
microlensing. To perform a complete analysis of this chromatic event our observations were used together 
with Optical Gravitational Lensing Experiment light curves. Both data sets cannot be reproduced by the simple 
phenomenology described under the caustic crossing approximation; using more realistic representations of 
microlensing at high optical depth, we found solutions consistent with simple thin disk models (r, a \^^^); 
however, other accretion disk size-wavelength relationships also lead to good solutions. New chromatic events 
from the ongoing narrow band photometric monitoring of Q 2237+0305 are needed to accurately constrain the 
physical properties of the accretion disk for this system. 

Subject headings: gravitational lensing — accretion, accretion disks — quasars: individual (Q 2237+0305) 



1. INTRODUCTION 

Gravitational lensing is independent of wavelength 
(Schneider, Ehlers & Falco 1992). However, in many 
gravitationally lensed quasars, differences in color between 
the images are observed. These chromatic variations could 
be produced by two effects: differential extinction in the lens 
galaxy and chromatic microlensing. When each image's light 
crosses the interstellar medium of the lens galaxy it may be 
affected in different amounts by patchily distributed dust. 
This results in differential extinction between pairs of images, 
and makes possible the determination of the extinction law 
of the lens galaxy (Nadeau et al. 1991; Falco et al. 1999; 
Motta et al. 2002; Munoz et al. 2004; Mediavilla et al. 2005; 
EKasdottir et al. 2006). The other chromatic phenomenon 
arises when stars or compact objects in the lens galaxy are 
nearly aligned with the line of sight between the quasar image 
and the observer. Due to the relative motion between the 
quasar, the lens and the observer, the quasar image undergoes 
a magnification or demagnification known as microlensing 
(Schneider, Kochanek & Wambsganss 2006 and references 
therein). Therefore, fluctuations in the brightness of a 
quasar image will be a combination of the intrinsic quasar 
variability and the microlensing from the stars or compact 
objects in the lens galaxy. The intrinsic variability of the 
lensed quasar will appear in all images with certain time 
delay due to the different light travel times. Once this delay 
is determined, the light curves of the different images can 
be shifted, and then subtracted. The remaining fluctuations 
can be assumed to be caused only by microlensing. The 
microlensing magnification depends on the angular size of 
the source, in this case on the accretion disk of the quasar. 
Because the accretion disk is hotter closer to the black hole, 
and because the emission of the accretion disk depends on 
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temperature, different magnifications may be observed at 
different wavelengths. This effect is known as chromatic 
microlensing, and it offers unprecedented perspectives into 
the physical properties of accretion disks (Wambsganss & 
Paczynski 1991). Its detection in a lens system will lead to 
accurate constraints in the size-wavelength scaling. 

Different authors have made some attempts to detect chro- 
matic microlensing and use it with different scopes. However, 
in many cases the detection of chromatic microlensing was 
rather ambiguous, since in general it is not easy to disentan- 
gle this effect from others with similar observational signa- 
tures (Wisotzki et al. 1993; Nadeau et al. 1999; Wucknitz 
et al. 2003; Nakos et al. 2005). The first significant appli- 
cations of chromatic microlensing appeared in a very recent 
work of Poindexter et al. (2008). They have used chromatic 
microlensing in order to determine the size of the accretion 
disk of HE 1 104-1805 as well as its size-wavelength scaling. 
More recently, Anguita et al. (2008) have also made these 
kinds of studies in Q 2237+0305. In both cases they found 
solutions consistent with the simple thin disk model (Shakura 
& Sunyaev 1973), but stronger constraints would be needed 
to accurately determine the size of the accretion disk. 

In this work, we present a chromatic microlensing detec- 
tion in one of the images of Q 2237+0305 (Huchra et al. 
1985). This lens system has very good properties to study 
microlensing events. The light coming from a distant quasar 
at = 1 .695 is deflected by a nearby spiral galaxy {zl = 0.039) 
forming four quasar images nearly symmetrically distributed 
around the lens center. Since the light passes through the cen- 
tral part of the galaxy, where the microlensing optical depth is 
very large, high magnification events (HMEs) occur very fre- 
quently. Moreover, since the galaxy has a very low redshift, 
and due to the symmetry in the positions of the images, time 
delays are expected to be very small (< 1 day; Vakulik et al. 
2006; Koptelova et al. 2006). Therefore, time delay correc- 
tions are not needed to study microlensing in this lens system, 
and the flux ratio between the images will remove the contri- 
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TABLE 1 
Log of ALSFOC Observations 





Observation Date 


Filter 


JjAUtJSLIlC \Jj 


Q 2237+0305 


2003 Aug 


26 


Str-u (A = 


= 3510 A) 


4x600 


Q 2237+0305 


2003 Aug 


26 


Str-v (A = 


=4110 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


Str-b (A = 


= 4670 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


Str-y (A = 


: 5470 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


Iac#28 (A 


= 6062 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


Ha (A = 


6567 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


Iac#29 (A 


= 7015 A) 


3 X 600 


Q 2237+0305 


2003 Aug 


26 


I-band (A 


= 8130 A) 


4x 100 


Q 2237+0305 


2003 Aug 28 


Str-u (A = 


= 3510 A) 


3 X 600 












2x900 


Q 2237+0305 


2003 Aug 


28 


Str-v (A = 


= 4110 A) 


3x600 


Q 2237+0305 


2003 Aug 


28 


Str-b (A = 


= 4670 A) 


3x420 


Q 2237+0305 


2003 Aug 


28 


Str-y (A = 


= 5470 A) 


3x420 


Q 2237+0305 


2003 Aug 


28 


Iac#28 (A 


= 6062 A) 


3x420 


Q 2237+0305 


2003 Aug 28 


Ha (A = 


6567 A) 


3x420 


Q 2237+0305 


2003 Aug 28 


Iac#29 (A 


= 7015 A) 


3x420 


Q 2237+0305 


2003 Aug 


28 


I-band (A 


= 8130 A) 


3 X 100 



2 - 



CIV 



Lya 



bution of the intrinsic quasar variability in the Hght curves. 

Specifically, we analyze a particular chromatic microlens- 
ing event, and study the physical scenarios in which this effect 
can be produced. The results obtained in this work are the first 
step toward obtaining precise constraints in the physical prop- 
erties of the quasar accretion disk. In §2 the observations and 
the data analysis are presented. The data fitting techniques 
and different models to describe the observed chromatic mi- 
crolensing are discussed in §3. A summary of the main con- 
clusions appears in §4. 

2. OBSERVATIONS AND DATA ANALYSIS 

On the nights of August 26 and 28 2003, we observed 
Q 2237H-0305 with the 2.56 m Nordic Optical Telescope 
(NOT) located at the Roque de los Muchachos Observatory, 
La Palma (Spain), using the 2048x2048 ALSFOC detector. 
Its spatial scale is 0.188 arsec/pixel. Seven narrow filters plus 
the wide I-Bessel were used. The whole set covered the wave- 
length interval 3510-8130 A. Table 1 shows a log of our ob- 
servations. 

Among the filters that we used only two were affected by 
the emission lines of the quasar. The Stromgren-u filter was 
affected by almost 40% by the Lya emission line, where the 
emission lines of the quasar and its nearby continuum were 
modeled according to the SDSS quasar composite spectrum 
(Vanden Berk et al. 2001) to perform this estimation. At 
A = 41 lOA the wavelength range covered by the Stromgren-v 
filter coincides with the position and width of the CIV emis- 
sion line. 

The data were reduced using standard procedures with 
IRAF packages, and PSF photometry fitting was used to de- 
rive the difference in magnitude versus wavelength between 
the four images of Q 2237-1-0305. The galaxy bulge was mod- 
eled with a de Vaucouleurs profile, and the quasar images 
as point sources. This model was convolved with different 
PSFs observed simultaneously with the lens system in each of 
the frames, and compared to the image through statistics ( 
McLeod et al. 1998; Lehar et al. 2000). Due to the good see- 
ing conditions (0".6 in I band), the results of the photometry 
were excellent even in the bluest filters. 

Table 2 presents the measurements of the differences in 
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Fig. 1. — Magnitude differences as a function of the observed wavelength 
for Q 2237+0305 (NOT data). Images A and C show a clear signature of 
microlensing in the continuum but not in the emission lines. Moreover, the 
chromatic variation observed in image A could only be explained as a con- 
sequence of chromatic microlensing. The wavelength widths covered by the 
emission lines that affect two of our filters (CIV and Lya) are indicated. 



TABLE 2 
Q 2237+0305 PHOTOMETRY 



Filter 


mg-niA 


mc-rriA 


mD-rriA 


Str-u (A = 3510 A) 


1.57±0.06 


1.57±0.06 


1.52±0.06 


Str-v (A = 41 10 A) 


1.41 ±0.02 


1.45±0.04 


1.39±0.02 


Str-b (A = 4670 A) 


1.67±0.03 


1.47±0.04 


1.59±0.04 


Str-y (A = 5470 A) 


1.65±0.02 


1.45±0.06 


1.58±0.03 


Iac#28 (A = 6062 A) 


1.58±0.03 


1.39±0.04 


1.52±0.04 


Ha (A = 6567 A) 


1.46±0.02 


1.32±0.03 


1.46±0.05 


Iac#29 (A = 7015 A) 


1.56±0.03 


1.38±0.02 


1.52±0.02 


I-band (A = 8130 A) 


1.40±0.03 


1.22±0.02 


1.38±0.04 



magnitude in each filter. In Figure [H we have plotted six 
magnitude differences versus wavelength (three of which are 
independent) between the four images of the quasar. 

The difference in magnitude between B and D does not 
show any chromatic variation, and it is reasonable to suppose 
that neither B nor D is undergoing flux variations with wave- 
length. The differences m^-mc, and mo-nic, show that both 
curves have the same trend, and this implies that the magni- 
tude of image C varies with wavelength (consistent with our 
hypothesis that neither image B nor D suffers of chromatic- 
ity). In these curves it is also observed that there is no magni- 
tude variation with the wavelength in the region between 4670 
A and 8130 A, but at the wavelengths corresponding to the 
Stromgren-v and Stromgren-u filters, which are the only ones 
affected by emission lines, we observed a deviation from the 
trend. As we stated above the Stromgren-u filter was affected 
by almost 40% by the Lya emission line, and at A = 4110 
A, the wavelength range covered by the Stromgren-v filter 
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Fig. 2. — Magnitude differences between Q 2237+0305 images obtained from OGLE V-band light curves. The dashed line indicates our observation date, and 
the solid lines limit the selected time interval in which the study was performed. 



coincides with the position and width of the CIV emission 
line. Then, the dependence of iriB-nic or mo-mc with wave- 
length observed under our filter set configuration is the ob- 
servational signature expected for image C being affected by 
microlensing in the continuum but not in the emission lines. 
The fact that microlensing is absent in the emission lines con- 
firms the relatively large size of the broad-line region (BLR) 
in Q 2237+0305 according to the existing relation between 
the size of the BLR and the intrinsic luminosity of the quasar 
(Kaspi et al. 2000, Abajas et al. 2002). 

Finally, we analyzed the differences between image A and 
the other ones. Again, as the trend of the curves is the same 
for each difference, we conclude that A is the image which 
undergoes flux variations. As the points corresponding to 
Stromgren-v and Stromgren-u filters deviate from the respec- 
tive trend, we can conclude again that image A is undergo- 
ing microlensing in the continuum but not in the emission 
lines. An even more interesting result appears from analyz- 
ing the points corresponding to the continuum, where image 
A shows a significant variation with wavelength. In principle, 
this effect could be produced by extinction and/or chromatic 
microlensing. But as image A is the one undergoing the flux 
variation, if it were due to extinction the slope of the curve 
should be negative, and not positive as it is observed in the 
plot. Therefore, we have strong evidence of chromatic mi- 
crolensing in image A. 



Since we have observations of only two consecutive nights 
in that period, to know what was the signature of the mi- 
crolensing in the system, we used the public Optical Grav- 
itational Lensing Experiment(OGLE) ^ data (Wozniak et al. 
2000). Since 1999 January OGLE has been monitoring 
Q 2237+0305 in the V-band. In Figure |2] three independent 
magnitude differences between the quasar images obtained 
from their data are plotted. In the time interval between HJD 
2452764 and HJD 2453145 (indicated in Figure |2] between 
solid lines) the differences (m^-mc) and {niB-mo) are almost 
constant. This would imply that the flux variations of the im- 
ages B, C, and D are caused mainly by intrinsic variability 
of the quasar, which agrees with our previous hypothesis that 
neither B, C nor D is undergoing chromaticity originated by 
an HME. The fact that in our observations we have detected 
(nonchromatic) microlensing in image C is also in agreement 
with OGLE data, since a constant micro-magnification value 
during this period justifies both observations. Therefore, the 
brightness fluctuations (~ 0.3 mag) observed in the difference 
(mg-m/i) can be assumed to be due only to microlensing- 
induced variability in image A. In other words, from Figure 
|2]we can say that image A was in fact undergoing an HME 
on the dates of our observations which correspond to HJD 
2452878.5 and 2452880.5. 

^ |http://bulge.princeton.edu/~ogIe7] 
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Fig. 3. — Typical observable effects of chromatic microlensing. The am- 
plitude of the microlensing is different at different wavelengths due to the 
expected dependence of the accretion disk size with wavelength. In addition, 
the maximum magnification is reached at a different instant with a time delay 
^tmax (see the text). 



Summarizing we can say that, based on the 2003 NOT ob- 
servations, images C and A are undergoing microlensing in 
the continuum but not in the emission Hnes, and that image A 
is also affected by chromatic microlensing. In the next sec- 
tion this unambiguous detection of chromatic microlensing is 
used, along with the OGLE data from HJD 2452764 to HJD 
2453145, to explore the phenomenology of the detected mi- 
crolensing event and to study the physical properties of the 
accretion disk. 

3. DATA FITTING 

3.1. Chromatic microlensing and caustic crossing 
approximation 

To review some basic aspects of chromatic microlensing 
we will start recalling the simplest model for microlensing, 
the caustic crossing approximation. Taking a Cartesian co- 
ordinate frame in which the caustic lies along the y-axis, the 
magnification at a distance x is proportional to the well-known 
expression \ + j3 H(x)/^ (Schneider & Weiss 1987). Here 
H represents the Heaviside step function, and (3 measures the 
caustic strength. Assume a brightness profile for the quasar 
accretion disk given by 



(1) 



where rj represents a typical size scale of the quasar accretion 
disk; provided a global factor (or a global constant in magni- 
tudes) containing the magnification background, the magnifi- 
cation of the quasar accretion disk is 



/3 



M~ 1 + ^0(0, 



(2) 



where ^ = ^ is the distance between the center of the disk 
and the caustic, xq, in units of the scale radius r^, and it can 
take negative or positive values to distinguish both sides of 
the caustic. Analytical expressions can be straightforwardly 
obtained for several accretion disk models, e.g., Gaussian 
(Schneider & Weiss 1987) or power laws (Shalyapin 2001). 
If a constant relative source plane velocity is considered be- 
tween the disk and the caustic, then jcq = v (f - fo), where fo is 
the instant of time at which the center of the disk is located 
at the caustic position. Therefore we can see that an HME 
due to a single caustic crossing is, in addition to the crossing 
time f(), a function of only two parameters: (3' = fij ^frl and 
r[ = Tj/v. Due to this degeneracy, the disk size and the ve- 
locity V cannot be determined independently but only the ratio 

Figure [3] shows the typical observable effects of chromatic 
microlensing. Because the size of the accretion disk is ex- 
pected to change with wavelength, the amplitudes of the HME 
events observed at different wavelengths are different. The 
maximum magnification, in magnitudes, relative to the local 
average microlensing magnification background (see Figure 
O will be 



Am„„.i. =-2.5 Log 



1 + 



(3) 



The value of 4){£,,mix) depends on the brightness profile model. 
For instance, for a Gaussian model 4){S,,„ax — 0.765) ~ 1.021. 
Hence, given a brightness profile model, the observable 
Am„,av provides a direct estimation of the dimensionless pa- 
rameter (3/ y/Fs. 

If an HME is monitored at two different wavelengths, then 
using equation (O the ratio between the different correspond- 
ing sizes is 

2 



10-'^"C,„/2-5_l 
„/2.5 _ I 



10- 



(4) 



and it is very remarkable that, provided that equation (|2|i is 
satisfied, this expression is independent of the brightness pro- 
file of the accretion disk. As an example we can use the re- 
cent results obtained by Anguita et al. (2008), where they 
analyzed another HME observed in the 1999 OGLE cam- 
paign. From Figure 2 in their paper it is possible to estimate 

that Am*^^. ~ -0.8 (see also OGLE El data for a better sam- 
pling) and Amjj,^^ ~ -0.57. Therefore, applying equation (|4|l 
r^/rf ^2.5. This value is slightly different from the ratio ob- 
tained in their paper where they use a magnification pattern 
instead of a caustic crossing. 

Another interesting property is obtained with the observ- 
able Atmax, i.e., the difference in time to reach the maxi- 
mum magnification for each source size (see Figure O. It 
is straightforward to obtain 



r^' 
r 



(5) 



Although it depends on the quasar disk model, it allows us 
to determine directly the size-velocity ratio. Unfortunately 
the light curve shown in Figure 2 of Anguita et al. (2008) 
has not the temporal sampling to accurately estimate Af,„n.r, 



where / is the modified Bessel function of the first kind. 
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Fig. 4. — OGLEV-band (me-mA) data. The best fit using a simple double- 
caustic crossing approximation and a Gaussian brightness profile is shown. 
The dashed line indicates our observation date. 

but assuming a value 10 days we would find the order of 
magnitude for r,* = r^/v ~ 20 day"' . 

3.2. Double -caustic crossing approximation 

We start using the simple model of the caustic crossing ap- 
proximation with the aim of finding a simultaneous fit of the 
difference {mB-niA) obtained from the OGLE V-band light 
curves (over the temporal range from HJD 2452764 to HJD 
2453145) (Figure nil and our chromatic observation. As ex- 
plained before, in this case, the fluctuations in (m^-m^) cor- 
respond only to microlensing-induced variability in image A. 
Because of the shape of the data we assumed that the image 
A of the quasar is undergoing a double-caustic crossing event, 
and we considered a Gaussian brightness profile^ for the ac- 
cretion disk I(r) = exp [-r^/2r^'\ . The main parameters 
involved in our fitting are: the ratio between r^ and the relative 
caustic-disk source plane velocity, r' = rs/v, the dimensionless 
caustic stretch divided by ^/r^, the time corresponding 

to each caustic crossing, f,=i.2, and the relative microlensing 
magnification background between each side of the caustics, 
where the subindex / = 1,2 refers to each of the caustics. 

The estimation of the model parameters was carried out 
with a minimization fitting method. Due to the scatter 
of the OGLE data, the points were binned to have only one 
point per day, which value corresponds to the averaged binned 
points, and which associated error corresponds to the disper- 
sion calculated in each day. Because the first caustic cross- 
ing event is not completely sampled, the parameters associ- 
ated with it cannot be determined; however the ones associ- 
ated with the second caustic can be bounded {(32 = 0.3 ±0.1, 
f2 = 3008 ± 12). Figure|4]shows the best fit that corresponds to 

^ Other brightness profiles were also used, as the thin accretion-disk ap- 
proximation (e.g. Kochanek 2004), or power laws (e.g. Shalyapin et al. 
2002), but the differences in the results produced by these other profiles are 
not significant. 



Fig. 5. — Image A "chromatic map". AMv represents the microlens- 
ing magnification in the V-band, and AMv-/ is the difference in magnitude 
mv -mi due to chromatic microlensing. The open square corresponds to the 
observed values. The different grayscales in the map coiTespond to different 
values of the probability of occurrence. 

/Nf = 0.22, where Nf is the number of degrees of freedom. 

Although this simple model produces a good light curve 
fitting for the OGLE data, it cannot reproduce the observed 
chromatic effect. Taking as reference the filters Stromgren-y* 
and I-band, we observe a chromaticity (m^-w^) = -0.24 ± 
0.05 mag . However, assuming a radius-wavelength scaling 
that goes as oc A'*/-' (thin-disk model), the chromaticity that 
can be produced by this double-caustic approximation model 
on the day of our observation (HJD 2452879.5) is - -0.01 
mag, and a maximum chromaticity of ^ -0.25 mag is only 
achieved for an unrealistic ratio rj/rj' > 10^. The limitation 
of the caustic crossing approximation in reproducing the ob- 
served chromatic microlensing magnification is due to the fact 
that the induced chromaticity is only a function of the caustic 
stretch j3 / ^JTs, and it is independent of the microamplification 
background. However, as we will see in the next subsection, 
this is not the case in a more realistic scenario. 

This is a good example which confirms, as stated by 
Kochanek (2004), that due to the high microlensing optical 
depth of Q 2237-1-0305 it is unlikely that the system under- 
goes a "clean" caustic crossing. In this sense the model that 
we used here is too simple to approach the real behavior of the 
system. The appropriate way to find a solution is to consider 
a more realistic physical scenario, i.e., a magnification pattern 
is needed to describe the observational signatures in this case. 

3.3. Modeling microlensing from magnification maps 

Magnification maps for images A and B would be needed 
to find trajectories that reproduce {ms - m^) V-band OGLE 
difference during the period of interest. As stated before, in 

* The central wavelengths of the Str-y and V-band filters almost coincide. 
Therefore hereafter we will consider Str-y and V-band filters as equivalents to 
simplify the comparisons between our data set and that obtained by OGLE. 
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Fig. 6. — Histogram of the number of tracks, Nj, that reproduce OGLE V- 
band data as well as the chromatic data set, as a function of the magnification 
ratio between images A and B, mg. 



the time interval between HJD 2452764 and HJD 2453145, 
the brightness fluctuations observed in (rtiB-mf^ are due only 
to microlensing in image A, since image B is mainly affected 
by intrinsic variability of the quasar. Therefore, the behav- 
ior of image B can be represented with a constant magnifica- 
tion value, and magnification patterns for this image are not 
needed. Its contribution to the (m^ - ni^ magnitude differ- 
ence, as well as any extinction correction or uncertainties in 
the macromodels will be included in a constant parameter, otq, 
that will shift image A microlensing fluctuations. 

We have built magnification patterns for image A using the 
inverse ray shooting technique (Wambsganss 1990, 1999) and 
the inverse polygon mapping (Mediavilla et al. 2006). It was 
assumed that all the mass responsible for microlensing is in 
compact objects of 1 M©, and the values for the convergence 
of K = 0.36 and for the shear of 7 = 0.4 (Schmidt et al. 1998) 
were considered. The map dimensions are of 4096 x 4096 
pixels which correspond to physical dimensions of 20 x 
20 ve- Therefore the magnification pattern has a resolution of 
0.005 r^/pixel. 

With the aim of evaluating the likelihood of reproducing 
the microlensing and chromatic microlensing magnifications 
we observed, a "chromatic map" for image A was built (Fig- 
ure Is). To each point in the map with coordinates (AMv, 
AMy-/), corresponds a value of the probability of occurrence. 
Here AMy = -2.5 log /xy represents the microlensing mag- 
nification in V-band, and AMy_/ is the difference in mag- 
nitude my- nil due to the disk size- wavelength dependence 
(i.e. chromatic microlensing). To compute those magnitude 
differences, since we are considering an extended object, the 
pattern was convolved with the brightness distribution of the 
source and, as before, we assumed a Gaussian profile. At 
A = 5430 A (V-band) a Gaussian width of r, = 0.03 (- 1 .5 
light days) was adopted which corresponds to the average 
value found by Kochanek (2004). At A = 8130 A (I-band), 



assuming that the accretion disk has a size-wavelength de- 
pendence that goes as r, oc A^/^ (thin-disk model), the pattern 
was convolved with a Gaussian width of rj = 0.05 3 
light days). 

Once the quantities AMy and AMy-i were evaluated at 
each pixel in the V-band convolved pattern and in the V-I dif- 
ference pattern respectively, the probability of occurrence was 
calculated, and its distribution is represented in Figure |5] 

As far as the baseline for no microlensing {i.e., the magnifi- 
cation ratio between the images in the absence of microlens- 
ing) is unknown; the true value of AMy is unknown. How- 
ever, our observations in the filters Str-u and Str-v which are 
affected by the emission lines would allow us, in principle, 
to estimate separately the true magnification ratio between 
the images, and the microlensing. This is because the ob- 
served flux for the image / is given by Fi = fi'""'^''"(ii1"'^''°Fc+Fi), 
where the subindex c refers to the continuum and / to the 
emission lines, i.e., the microlensing is affecting the contin- 
uum but not the emission lines. In practice this calculation 
is complicated to do because a model for the continuum and 
the emission lines have to be assumed, and also because small 
variations in the measured magnitudes introduce large uncer- 
tainties in both factors [p!^""'" and fi""^''"). For instance, these 
large uncertainties do not allow the determination of the mi- 
crolensing in image C. For image A, using our measurement 
in the Str-u filter and modeling the Lya emission line and its 
nearby continuum according to the SDSS quasar composite 
spectrum (Vanden Berk et al. 2001), we estimate a flux ratio 
Fb/Fa ^ Fd/Fa^ 0.6 ±0.2 and a microlensing magnification 
of /iy ~ -1.2 ± 0.4 mag. These flux ratios are in agreement 
with those found by Wayth et al. (2005). Moreover a similar 
macroamplification for the four images would be a reasonable 
scenario for a macrolens model describing the large symme- 
try observed in the four image positions, but it would imply 
that the anomalous brightness of the image A is produced by 
a large microamplification background. 

We found that given a microlensing magnification of = 
-1.2 ±0.4 mag for image A, the probability of observing a 
chromatic microlensing magnification m^-m'j^ = -0.24±0.05 
mag is of ^ 4%. Even though this value is not very high, 
the probability of finding a solution is now different from 
zero. That is, the modeling of microlensing chromaticity from 
microlensing magnification maps leads to a phenomenology 
(richer than that associated with the caustic crossing approx- 
imation) that includes the observed data. The square in the 
"chromatic map"(Figure|5]l corresponds to these values. The 
probability increases if, instead of using the cc X^/^ law, 
greater values of /r^ are considered, reaching a maximum 
value of ~ 27% for rj/rf ~ 3.9. 

The microamplification estimation of /iy = -1 .2 ± 0.4 mag 
could be affected by large uncertainties; for instance a real 
source spectrum with deviations from the Gaussian approxi- 
mation of the emission line profile (Vanden Berk al. 2001), 
or dust in the lens galaxy, would lead to different flux ratios 
between the images and microamplification values, and there- 
fore it must be only considered as a tentative estimation. A 
spectroscopic monitoring of the system would be needed to 
clarify the true magnification ratios and microlensing for the 
images. 

To complete our study we explore the possibilities of jointly 
reproducing the mg - ma V-band OGLE data and our obser- 
vations. In the first place, we fit the OGLE data using the 
V-band convolved pattern, where 2 x 10'' random trajecto- 
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Fig. 7. — In the top panel, the best OGLE hght curve fitting (dashed line; x~ = 1") is shown together with the fitting that also reproduces the chromaticity 
(X^ = 13; solid line). In the lower panel the chromatic effects produced along each trajectory are compared, where the filled square represents the observed 
chromaticity. 




Fig. 8. — Left: Trajectory that best fits the OGLE data, but which does not reproduce the observed chromaticity. Right: Trajectory that best reproduces both 
data sets. Both sections of the magnification pattern have a physical size of 2.9 r£ X 2.9 r£. 



ries were traced, and compared with the OGLE difference 
niB-niA through a statistics. The source plane velocity, 
and the macromagnified image difference were varied in the 
ranges 0.0015 HJO"! < v < 0.0065 HJO"', and -1.08 
mag < niQ < 2.59 mag respectively. The wide range of the 
covered source plane velocities includes already existing lim- 
its (Kochanek 2004), and the bounds set for the macromag- 
nified image difference were chosen according to covering a 
large range of microlensing magnifications (from -2.48 mag 
to 1.19 mag ). From a total of 7.16 x 10'" simulated trajecto- 
ries we have found 2711,493 tracks with /Nf < 1 (where 



Nf = 57) which satisfactorily fit the OGLE data. 

The next step is to determine in how many of those tracks 



the observed chromaticity (m^ ~h^a - -0.245 ±0.05) is re- 
produced. We found that 6443 of the tracks which reproduce 
OGLE data also reproduce, at the I - a level, our chromatic 
observation. The probability distribution of mo considering 
those trajectories is bimodal (Figure |6]l, with peaks located 
at mo ^ 0.1 and mo ^ 0.9, which corresponds to microlens- 
ing magnification values of /iy ^ -1 -4 and fiy ^ -0.6 respec- 
tively. Solutions are found in all the velocity range selected 
above with the number of tracks increasing when the veloc- 
ity decreases. The best fit is obtained with vq = 0.003 
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Fig. 9. — Magnitude difference between the images A and B of 
Q 2237+030 as a function of wavelength. The filled squares represent the ob- 
served values that are not affected by emission hues, while the empty squares 
are associated with the two filters contaminated by emission lines. The soHd 
line connects the values obtained for the date of our observation considering 
the track that best fits both data sets. These values were calculated by con- 
volving image A pattern with different source sizes coiresponding to each 
of the filters we used, and we assumed an accretion disk size-wavelength 
scaling that goes as Ts ~ A**/' . For comparison, we plot (dashed line) a lo- 
cal fit to the data modeling the chromatic microamplification as a power law 
(/.(A) ~ A-«'^). 

HJD~' and niQ = 0.87 mag which impHes fiy = -0.62 mag. 
In Figure [T] (top panel) the OGLE fit using these parameters 
is shown (soHd Une; = 13), and compared with the best 
OGLE data fitting which does not reproduce the chromaticity 
(dashed Une; = 10)- The chromatic behavior along each 
of these tracks is compared in Figure |7] (bottom panel). As 
shown in Figure[8] both trajectories are located in a complex 
caustic configuration. 

Therefore we have found models that describe both data 
sets under a thin disk model profile. Finally to compensate 
for our lack of a temporal sampling we use the good resolu- 
tion in wavelength of our observations to check the reliability 
of the thin-disk model. Considering that the size-wavelength 
scaling goes as r, cc we have convolved the image A 
magnification pattern with Gaussian profiles of different r, 
corresponding to each of our filters. For the track that best 
fits the data sets, the chromatic microlensing magnification 
was computed for the day of our observation in each filter 
The values predicted by this model are in perfect agreement 
with the observed ones that correspond to filters not affected 
by emission lines (see Figure |9] solid line). This does not 
prove that the r, a X^l^ law is the best model fitting the data 
however, it is remarkable, that such a model is fully consistent 
with our observations. 

Needless to say, other accretion disk size-wavelength rela- 
tionships will also lead to good solutions of both data sets. In 



fact, bigger r[ / rj' ratios produce a larger number of tracks that 
describe the observed chromaticity and microlensing magni- 
fications correctly. Therefore we show that our observations 
are compatible with a thin-disk model, although laws with 
larger r{ / ratios would have higher probability because they 
require small microlensing backgrounds. In any case a sin- 
gle chromatic event does not allow us to constrain enough 
the variation of the source size with wavelength, what can 
only be done with several chromatic microlensing detections 
at different epochs, and performing a global fit to the data. 
For this reason we are developing an ongoing project where 
Q 2237-1-0305 is monitored with the same filter set used in this 
work. 



4. CONCLUSIONS 

Using excellent quality ground-based data (narrowband 
photometry taken with the NOT at a single epoch) we have de- 
tected microlensing in the continuum but not in the emission 
lines in two images of Q 2237-1-0305 (A and C). This detection 
confirms that in luminous quasars such as Q 2237-1-0305 the 
BLR is large enough to escape global microlensing magnifica- 
tion. The data unambiguously indicate that the microlensing 
undergoing by image A is chromatic. This is an interesting 
detection, since this phenomenology is a fundamental tool to 
understand the structure of quasar accretion disks. 

Using our data set together with the OGLE data, we have 
explored different physical scenarios in which both data sets 
can be fitted. The observations cannot be reproduced by the 
simple phenomenology described under the caustic crossing 
approximation. Modeling microlensing from magnification 
patterns and considering a Gaussian brightness profile for the 
accretion disk we found solutions compatible with a simple 
thin-disk model. The best solution matches remarkably well 
the variation of microlensing magnification with wavelength 
according to the r, oc law. 

However, the low probability of the solutions leads us to 
think that the considered microlensing/quasar models could 
not be enough complete or realistic to describe the phe- 
nomenology of chromaticity. In fact we have seen that the 
probability of the solutions increases for laws with larger 
rj/rj' ratios than the 4/3 law. The single chromatic event 
detected in Q 2237-1-0305 looks insufficient to explore more 
complex microlensing scenarios. Additional chromatic detec- 
tions (expected from our ongoing narrowband monitoring of 
Q 2237-1-0305) are needed to obtain an accurate constraint in 
the size-wavelength relation of the accretion disk for this lens 
system. 
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